Dynamical polarizability of atoms in arbitrary light fields: general theory and 

application to cesium 
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We present a systematic derivation of the dynamical polarizability and the ac Stark shift of 
the ground and excited states of atoms interacting with a far-off-resonance light field of arbitrary 
polarization. We calculate the scalar, vector, and tensor polarizabilities of atomic cesium using 
resonance wavelengths and reduced matrix elements for a large number of transitions. We analyze 
the properties of the fictitious magnetic field produced by the vector polarizability in conjunction 
with the ellipticity of the polarization of the light field. 
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I. INTRODUCTION 

One of the main motivations of current laser cooling 
and trapping techniques is to use atoms for storing and 
processing quantum information that is encoded in the 
atomic states by means of resonant or near-resonant light. 
Due to the weak coupling of neutral atoms to their envi- 
ronment, coherent manipulation of atomic states can be 
robust against external perturbations This makes 
optically trapped neutral atoms prime candidates for, 
e.g., the implementation of quantum memories and quan- 
tum repeaters @~iH' For atom trapping, far-off- resonance 
laser fields are used because they ensure low scattering 
rates, compatible with long coherence times. The pres- 
ence of these intense far-detuned light fields shifts the 
energy levels of the atom. In general, the light shift (ac 
Stark shift) depends not only on the dynamical polariz- 
ability of the atomic state and on the light intensity but 
also on the polarization of the field. For this reason, var- 
ious experimental situations require a systematic study 
of the dynamical polarizability of the ground and ex- 
cited states of atoms interacting with a far-off-resonance 
light field of arbitrary polarization. In particular, this 
becomes important for optical trapping using near-fields 
or nonparaxial light beams. One example is nanofiber- 
based atom traps, which have recently been realized [E|, [f| 
and in which the nanofiber-guided trapping light fields 
are evanescent waves in the fiber transverse plane 
Another example is tightly focused optical dipole traps, 
where the longitudinal polarization component of a non- 
paraxial light beam can lead to significant internal-state 
decoherence @, Q . Plasmonically enhanced optical fields 
PUSH also have, in general, complex local polarizations. 
Therefore, the calculation of the resulting optical poten- 
tials in all these cases requires a suitable formalism to 
take polarization effects into account. 

Despite a large number of works on the polarizabili- 
ties of atoms, most of the previous calculations were de- 



voted to the static limit [12-16]. Accurate polarizabil- 
ities for a number of atoms of the periodic table have 
been calculated by a variety of techniques [Hj]. These 
include the sum-over-states method, which is based on 
the use of available experimental and/or theoretical data, 
and the direct methods, which are based on ab ini- 
tio calculations of atomic wave functions. The ab ini- 
tio calculations of atomic structures involve the refined 
many-body perturbation theory, the relativistic coupled- 
cluster calculations, or the random phase method [l6| . 
High-precision ab initio calculations of atomic polariz- 
ability have been performed using the relativistic all- 
order method in which all single, double, and partial 
triple excitations of the Dirac-Fock wave functions are 
included to all orders of perturbation theory [3 fl5l - 
Recently, in order to search for magic wavelengths [17| 
for a far-off-resonance trap, the dynamical scalar and 
tensor polarizabilities as well as the light shifts of the 
ground and excited states of strontium [13, HH and ce- 
sium [l|| [2(| have been calculated for a wide range of 
light wavelengths. The principal idea of magic wave- 
lengths is based on a clever choice of the trapping light 
wavelength for which the excited and ground states of 
an atom experience shifts of equal sign and magnitude 
171 . Magic wavelengths have been found for atomic ce- 
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sium in red-detuned traps (l7l [l9| and in combined two- 
color (red- and blue-detuned) traps [20[. Searches for 
magic and tune-out wavelengths of a number of alkali- 
metal atoms (from Na to Cs) have been conducted by 
calculating dynamical polarizabilities using a relativistic 
coupled-cluster method [U [22j. All the three compo- 
nents of the dynamical polarizability, that is, the scalar, 
vector, and tensor polarizabilities [231 ] . and the associ- 
ated ac Stark shifts have been calculated for the cesium 
clock states [2~3 - [26| . Calculations of the adiabatic po- 
tentials for atomic cesium in far-off-resonance nanofiber- 
based traps [El [27i - [29l | have been performed [El. I20I . [28l - [30l | . 
The vector polarizability was omitted in [2oTl28lj . but was 
included in the calculations for the ac Stark shifts in Ref. 
[30j |. The scalar, vector, and tensor polarizabilities of 
atomic rubidium have recently been calculated [3l| . 
Due to the complexity of the calculations for the dy- 
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namical polarizability of a realistic multilevel atom, vari- 
ous approximations have been used and different expres- 
sions for the components of the dynamical polarizability 
have been presented in different treatments. One exam- 
ple is that the counter-rotating terms in the atom-field 
interaction Hamiltonian was neglected in Refs. [Ill [2(| 
but was taken into account in Refs. (20l - [23 |. Another 
example is that the definition for the reduced matrix 
element used in Refs. [20l - [24j is different from that in 
Refs. [2f| (2(| |3(|. Furthermore, the coupling between 
different hyperfine-structure (Ms) levels of the same fine- 
structure state was taken into account in Refs. [2(| Hl| 
but was neglected in Refs. [UH^I- In addition, the nu- 
merical calculations require the use of resonance wave- 
lengths and reduced matrix elements of a large number 
of atomic transitions, which are not available in a single 
source. Since the authors of previous works often did not 
describe in detail the formalisms and the data they used, 
it is not easy to see the connections between their results 
and to employ them correctly. 

The purpose of this article is to provide a systematic 
treatment of the dynamical polarizability of the ground 
and excited states of atoms interacting with a far-off- 
resonance light field of arbitrary polarization. We specify 
all theoretical definitions and tools necessary for com- 
puting the light shifts of atomic levels. Based on the 
approach of Rosenbusch et al. [Hj], we provide the de- 
tails of the derivation of the expressions for the ac Stark 
interaction operator and for the scalar, vector, and ten- 
sor components of the dynamical polarizability. We also 
discuss the light-induced fictitious magnetic field. We 
supply a comprehensive set of experimental and theoret- 
ical data for resonance wavelengths and reduced matrix 
elements for a large number of atomic transitions that 
allows one to perform the computation of the light shifts 
of the levels associated with the r>2-line transition of ce- 
sium. Furthermore, we present the results of numerical 
calculations for the corresponding components of the po- 
larizability for a wide range of light wavelengths. Both, 
the atomic data and the numerical results are provided 
as electronic files which accompany this article [32J. 



A. General theory 

1. Hyperfine interaction 

We consider a multilevel atom. We use an arbitrary 
Cartesian coordinate frame {x,y,z}, with z being the 
quantization axis. In this coordinate frame, we specify 
bare basis states of the atom (see Fig. [T] for the levels 
associated with the D2-lhie transition of cesium). Due to 
the hfs interaction, the total electronic angular momen- 
tum J is coupled to the nuclear spin I. The hfs interaction 
is described by the operator [l[ 
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hA his I • J + frB] 
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Here, Ahf s and Bhfs are the hfs constants. Note that Ahf s 
and i?hfs depend on the fine-structure level \nJ). In the 
case of atomic cesium, the values of these constants are 
A hfs /2ir = 2298.1579425 MHz H and B hfs /2n = for 
the ground state 6S 1/2 and A his /2n = 50.28827 MHz and 
B his /2n = -0.4934 MHz [H for the excited state 6P 3/2 . 
We also note that high-order hfs interaction effects, which 
mix different fine-structure levels |nJ), have been omitted 
in expression ([1]) for the hfs interaction operator V hfs . 

Due to the hfs interaction, the projection J z of the to- 
tal electronic angular momentum J onto the quantization 
axis z is not conserved. However, in the absence of the 
external light field, the projection F z of the total angu- 
lar momentum of the atom, described by the operator 
F = J + 1, onto the quantization axis z is conserved. We 
use the notation \nJFM) for the atomic hfs basis (F ba- 
sis) states, where F is the quantum number for the total 
angular momentum F of the atom, M is the quantum 
number for the projection F z of F onto the quantization 
axis z, J is the quantum number for the total angular 
momentum J of the electron, and n is the set of the re- 
maining quantum numbers {nLSI}, with L and S being 
the quantum numbers for the total orbital angular mo- 
mentum and the total spin of the electrons, respectively. 
In the hfs basis {\nJFM)}, the operator I^ hfs is diagonal. 
The nonzero matrix elements of this operator are 

1 
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G(G + 1) - 2I{I+ 1)J(J+ 1) 



II. AC STARK SHIFT AND ATOMIC 
POLARIZABILITY 



21(21 - 1)2 J(2 J- 1) 
where G = F(F + 1) - 1(1 + 1) - J(J + 1). 



(2) 



In this section, we present the basic expressions for the 
ac Stark shift operator and the scalar, vector, and tensor 
polarizabilities of a multilevel atom interacting with a 
far-off-resonance light field of arbitrary polarization [23l — 
[26] ] . We also provide the results of numerical calculations 
for atomic cesium for a wide range of light wavelengths. 



2. ac Stark interaction 

Consider the interaction of the atom with a classical 
light field 

E = ]-Se- iUt + c.c. = \£ue~ iut + c.c, (3) 
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UJ b - UJ a - U - Hba/2 
1 

UJb - LO a + <JJ + i^ba/Z 



\b)(b\, 
\b)(b\. (7) 



F=3 



M= -3 



M=3 



We assume that V EE is the operator for the ac Stark in- 
teraction dHHH, i.e., that it correctly describes not only 
the level shift but also the level mixing of nondegenerate 
as well as degenerate states. While this educated guess 
has not been derived from first principles, it is consistent 
with the results of the second-order perturbation theory 
for the dc Stark shift [12, Il3ll and of the Floquet formal- 
ism for the ac Stark shift [23, [U| . 



FIG. 1. Energy levels associated with the D2 line of a cesium 
atom. 



where uj is the angular frequency and £ = £u is the 
positive-frequency electric field envelope, with £ and u 
being the field amplitude and the polarization vector, 
respectively. In general, £ is a complex scalar and u is a 
complex unit vector. 

We assume that the light field is far from resonance 
with the atom. In addition, we assume that J is a good 
quantum number. This means that we treat only the 
cases where the Stark interaction energy is small com- 
pared to the fine structure splitting. In the dipole ap- 
proximation, the interaction between the light field and 
the atom can be described by the operator 



V E = E d = -£ude 



-£*u* • de 1 "*, 
2 



(4) 



where d is the operator for the electric dipole of the atom. 
When the light field is far from resonance with the atom, 
the second-order ac Stark shift of a nondegenerate atomic 
energy level \a) is, as shown in Appendix \K\ given by 



5E a 



fee 



lu-d|a)l 2 



UJ b - UJ a 



ilba/2 



(a|u-d|6)p 



UJ b - UJ a + LU + i^ ba /2 



(5) 



Here, |a) and \b) are the atomic eigenstates with un- 
perturbed energies Tiui a and Tuob, respectively, and with 
spontaneous decay rates "f a and 7&, respectively, while 
Iba = 7a + 76 is the transition linewidth. We can 
consider the energy shift ([5]) as an expectation value 
8E a = (a\V EE \a), where 



V 



EE 



\£f 



[(u* ■ d)ft+(u • d) + (u ■ d)7Mu* • d)], (6) 



3. Atomic polarizability 

Let us examine the energy shifts of levels of a single 
fine-structure state \nJ). In general, due to the degener- 
acy of atomic levels and the possibility of level mixing, 
we must diagonalize the interaction Hamiltonian in order 
to find the energy level shifts. Since the atomic energy 
levels are perturbed by the Stark interaction and the hfs 
interaction, the combined interaction Hamiltonian is 



V 



hfs 



V 



EE 



(8) 



In terms of the hfs basis states | (nJ)FM) = \nJFM), the 
Stark operator V EE , given by Eq. ©, can be written as 



yEE = 



E 

FMF'M' 



V^ IF , M ' I (nJ)FM)((nJ)F'M'\, (9) 



where V^ IF , M , = {{nJ)FM\V EE \inJ)F l ' M') are the 
matrix elements and are given as [2J| 



V EE 



K = 0,1,2 
q=-K t . .. ,K 



"ij^ 11 * ® U }Kq 



x (_i)WJf+ j -J/ v /( 2 f + i)(2F' + 1) 

(10) 
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Here we have introduced the notations 
«nJ = {-1) K+W V2K + 1 

E(-D r { J r ) }l("' J '» d l 



nJ)f 



1 



x -Rel 

n \Un'J'nJ - W - iln'J'nj/2 



(-1) 



A" 



Un'J'nJ + W + i"fn'J'nj/2 



(11) 



with K = 0,1, 2, for the reduced dynamical scalar (K = 
0), vector (K = 1), and tensor (K = 2) polarizabilities 
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of the atom in the fine-structure level \nJ). In Eqs. (|T0|) 
and (fTTj) , we have employed the notations ( £ ^ 2 ^ ) 
and { 3 } j 3 6 } for the Wigner 3-j and 6-j symbols, re- 
spectively. The notations u n iji n j = 0J n 'ji — uj n j and 
Jn'J'nJ = In' j' + Jnj stand for the angular frequency 
and linewidth, respectively, of the transition between the 
fine-structure levels \n'J') and \nJ). The details of the 
derivation of Eqs. (|10[) and (fTT|) are given in Appendix 
151 Note that the above-defined polarizabilities are just 
the real parts of the complex polarizabilities. The imag- 
inary parts of the complex polarizabilities are related to 
the scattering rate of the atom (36j | . 

The compound tensor components {u* £g) u}xq in Eq. 
(ITOl) are defined as 



{U* ® u} K q = 



(-1)*+" 



E 

iu,/i'=0,±l 

xV2F+lf 1 K \\ (12) 

Here, u_i = (u x — iu v )/y/2, uo = u z , and ui = —{u x + 
iu y )/\/2 are the spherical tensor components of the po- 
larization vector u in the Cartesian coordinate frame 
{x,y,z}. 

The reduced matrix elements (n' J'\\d\\nJ) of the elec- 
tric dipole in Eq. (TTT|) can be obtained from the oscillator 
strengths 



fnjn' J' 



2m e w n > 



e^n' J'nJ 



1 



3?ie 2 2J+1 



|(n'J'||d||nJ)| 2 , (13) 



where m e is the mass of the electron and e is the elemen- 
tary charge, or from the transition probability coefficients 



A n , v 



n' J'nJ 



^n'J'nJ 



1 



3TTe hc 3 23' + 1 



|(n'J'||d||nJ)| 5 



(14) 



We note that the Stark interaction operator © with 
the matrix elements (ITU)) can be written in the form (23l . 



V EE = -\\£\ 



T 3[(u* • J)(u • J) + (u • J)(u* • J)] - 2J 2 



. v [u* x u] • J 



l nJ- 



23(23- 1) 



.(15) 



Here, a* j, cx^jj and a^j are the conventional dynamical 
scalar, vector, and tensor polarizabilities, respectively, of 
the atom in the fine-structure level |nJ). They are given 
as 24] 



1 



a 



nj 



l nJ 



* (0) 

l nJi 



23 



(J + 1)(2J+1) 



a nJ = -\ 



23(23-1) 



3(3+ 1)(2J+ l)(2J + 3) 



l nJ- 



(16) 



Note that for J = 1/2 and K = 2, the Wigner 6-j 
symbol in Eq. (fTI]) is zero. Thus, the tensor polarizabil- 
ity vanishes for 3=1/2 states (e.g., the ground states 
of alkali-metal atoms). In the case of linearly polarized 
light, the polarization vector u can be taken as a real 
vector. In this case, the vector product [u* x u] vanishes, 
making the contribution of the vector polarizability to 
the ac Stark shift to be zero. We also note that j n 'j'nj 
can be omitted from the denominators in Eqs. ([5]), (J7J, 
and pip when the light field frequency uj is far from res- 
onance with the atomic transition frequencies uj n 'j' n j. 

In general, V EE is not diagonal neither in F and nor 
in M. Therefore, in order to find the new eigenstates and 
eigenvalues, one has to diagonalize the Hamiltonian ©, 
which includes both the hfs splitting and the ac Stark 
interaction. However, in the case where the Stark inter- 
action energy is small compared to the hfs splitting, we 
can neglect the mixing of atomic energy levels with differ- 
ent quantum numbers F. In this case, the Stark operator 
yEE f Qr -j-jjg a t om m a particular hfs level \n3F) can be 
presented in the form [261 ] 



V EE = --\£\ z {a 



1 



2 J „s 



41-1 | ^nJF la nJF c^p 

j, 3[(u* • F)(u ■ F) + (u • F)(u* ■ F)] — 2F 2 



[u* x u] ■ F 



+ a nJF 

where 



a 



nJF — a nJ 



2F(2F-1) 



x nJ' 



. (17) 



V3(2J + 1) 



l nJF 



l nJF 



(-iy 



J+I+F 2F(2F + 1)(F I F \ a q) 



F + l 



3i3 



-(-1) 



J+I+F 



F 2 F 
3 13 



t 2F(2F- 1)(2F+ 1) 
3(F+ l)(2F + 3) 



(18) 



The coefficients a s nJFl a 



nJF and a^jF are the conven- 
tional scalar, vector, and tensor polarizabilities of the 
atom, respectively, in a particular hfs level. Note that 
the scalar polarizability ot s nJF does not depend on F. 
This statement holds true only in the framework of our 
formalism, where the hfs splitting is omitted in the ex- 
pression for the atomic transition frequency uv j> F'njF 
in the calculations for the atomic polarizability, that is, 
where the approximation 0J n 'j'F'njF = ^n'J'nJ is used. 
We also note that, if energies including hfs splittings are 
used in the denominators in the perturbation expression 
([5]), then the wave functions of the states \a) and \b) in 
the numerators should also incorporate hfs corrections to 
all orders of perturbation theory [24j, |37| . We emphasize 
that Eq. (fT7|) is valid only when the coupling between 
different hfs levels \n3F) is negligible. Thus, Eq. (fP7j) is 
less rigorous than Eq. (fT5"]) . 

Furthermore, we note that, when the off-diagonal cou- 
pling is much smaller than the Zeeman splittings pro- 
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duced by an external magnetic field B , the mixing of dif- 
ferent Zeeman sublevels can be discarded. In this case, 
the ac Stark shift of a Zeeman sublevel \FM) (specified 
in the quantization coordinate frame {x, y, z} with the 
axis z parallel to the direction zb of the magnetic field 
B) is given by 



AK 



EE 

FMFM 



Da T 



1 

3M 2 



l n.JF 



Ca 



nJF 



M 
2F 



F(F + 1) 



nJF~ 



2F(2F - 1) 



where 



C = |u_i| 2 — |ui| 2 = 2Im(u*M w ), 
D = 1-3KI 2 = l-3|u z | 2 . 



(19) 



(20) 



The coefficients C and D are determined by the polar- 
ization vector u of the light field at the position of the 
atom. Note that the parameter C, which characterizes 
the vector Stark shifts, depends on the ellipticity of the 
light field in the transverse plane (x,y). This parameter 
achieves its maximal magnitude |C| = 1 when the lon- 
gitudinal component of the field is absent and the light 
field is circularly polarized in the plane (x,y). We also 
note that the parameter D, which characterizes the ten- 
sor Stark shifts, vanishes when \u z \ ~ l/\/3. 



4- Fictitious magnetic field 

It is clear from Eqs. (|T5|) and (fl~7|) that the effect of 
the vector polarizability on the Stark shift is equivalent 
to that of a magnetic field with the induction vector jsF 
H 



B 



Act 



IS* x£} = 



l n.JF 



SflBgnJFF 



i[£*x£}. (21) 



Here, /j,b is the Bohr magneton and g n j and g n jF are the 
Lande factors for the fine-structure level \nJ) and the hfs 
level \nJF), respectively. The nonrelativistic value of the 
Lande factor g n j is given by [l[ 



g n j — gh - 



J(J + l)+L(L + l)-S(S+l) 



gs- 



2J(J + 1) 
J(J+1) + S(S+1)-L(L+1) 
2J(J + 1) 



(22) 



Here, gt = 1 and gs — 2.0023193 are the orbital and spin 
g-factors for the electron, respectively. When the contri- 
bution of the nuclear magnetic moment is neglected, the 
Lande factor g n jF is 



g n jF = g n j- 



F(F+1) + J(J + 1) - 7(1 + 1) 



2F(F + l) 



(23) 



The direction of the light-induced fictitious magnetic 
field B fict is determined by the vector i[£* x £], which is a 
real vector. Similar to a real magnetic field, the fictitious 



magnetic field B fict is a pseudovector, that is, B fict does 
not flip under space reflection. Another similarity is that 
both the real and fictitious magnetic fields flip under time 
reversal. If the light field is linearly polarized, we have 
i[£* x £} = and hence B = 0. The middle expression 
in Eq. (f2Tj) shows that B fict is independent of F, that 
is, B fict is the same for all hfs levels \nJF) of a fine- 
structure level \nJ). Comparison between the middle 
and last expressions in Eq. (12T1) shows that the factor 
OL v nJF l g n jpF does not depend on F. This conclusion is 
consistent with the relation 



l nJF 



F(F + 1) + J(J + 
(F + 1)^2J(J + 1)(2J + 1) 



(24) 



which can be obtained directly from the second expres- 
sion in Eqs. (|18[) with the use of an explicit expression 



for the Wigner 6-j symbol { j } F r } . 

In general, the vector Stark shift operator can be ex- 
pressed in terms of the operator J as 



EE 



fict^ 



(25) 



In the special case where the mixing of different hfs levels 
is negligible, that is, when F is a good quantum number, 
the vector Stark shift operator can be expressed in terms 
of the operator F as 



EE 



H B g n jF(F-B Rct ) 



(26) 



The vector form of Eqs. (j2l>]) and (f2"6")l allows us to con- 
clude that the fictitious magnetic field B fict can be sim- 
ply added to a real static magnetic field B if the latter 
is present in the system [see Eqs. (|C1[) and (|C3[) in Ap- 
pendix [C] . 

Let us discuss the case of the ground state nS\/2 of an 
alkali- metal atom. In this case, we have J = 1/2 and, 
therefore, a^j = 0. We assume that the hfs splitting 
of the ground state is very large compared to the Stark 
interaction energy. Then, the mixing of two different hfs 
levels _F = /±l/2of the ground state can be neglected, 
that is, F can be considered as a good quantum number. 
It is obvious that M is also a good quantum number when 
the quantization axis z coincides with the direction of the 
fictitious magnetic field B fict . 

For the hfs levels F = 7+1/2 of the ground state n5i/ 2 , 
we have g n jF\F=i+i/2 = -flW-pl^w-i/a = 5«j/(2I + 
1). When the hfs splitting of the ground state is very 
large compared to the light shift, the vector Stark shift 
operator is given in terms of the operator F by Eq. (|26[) . 
Hence, when the direction of the fictitious magnetic field 
B fict is taken as the quantization axis z, the vector Stark 
shifts of the sublevels M of the hfs levels F = I + 1/2 
and F = I — 1/2 of the ground state are 



EE i 



\F=I+l/2 



V-BgnJ 
21+1 



MB 



fict 



and 



EE I 



\F=I-l/2 



"21 + 1 



(27) 



(28) 
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respectively. These shifts are integer multiples of the 
quantity /i_B5njP fict /{2I + 1). In other words, as ex- 
pected from analogy with the well-known Zeeman effect, 
the shifts are equidistant with respect to the quantum 
number M . It is clear that the sublevels M and —M of 
the hfs levels F = 1+ 1/2 and F = I — 1/2, respectively, 
of the ground state have the same vector Stark shift. In 
contrast, the sublevels with the same number M of two 
different hfs levels F = I ±1/2 have opposite vector Stark 
shifts. Since the scalar Stark shift does not depend on F, 
the differential shift of the energies of the sublevels M' 
and M of the hfs levels F' = I + 1/2 and F = I - 1/2, 
respectively, of the ground state is just the differential 
vector Stark shift and is given by 

AWW=f^(M' + M)B fict 

= 8J§TT) {M ' + MMS * X£]1 (29) 

This differential shift vanishes when M' + M = 0. This 
result is valid only in the framework of our formalism 
where the hfs splitting is neglected in the calculations for 
the atomic polarizability. 



B. Numerical calculations 

We now present the results of numerical calculations 
for the dynamical scalar, vector, and tensor polarizabil- 
ities of the ground and excited states associated with 
the P/2-line transition of atomic cesium. Before we pro- 
ceed, we note that, in order to search for red- and blue- 
detuned magic wavelengths for a far-off-resonance trap, 
the scalar and tensor polarizabilities of the ground and 
excited states of atomic cesium have been calculated fl9i — 
l2lj . Releva nt p arameters were taken from a number of 
sources [lj| l49l - l5l| . Very recently, the vector light shifts 
of cesium atoms in a nanofiber-based trap have been 
studied [3(| • However, the results for the vector polariz- 
ability have not been explicitly provided. 

Our calculations for the polarizabilities of cesium are 
based on Eqs. (fTTH) in conjunction with Eqs. (fTTj) . The 
calculations for the polarizability of the ground state 
6S , 1 / 2 incorporate the couplings 65 , 1 / 2 (6-40)Pi/ 2i3 / 2 . 
The calculations for the polarizability of the excited state 
6P 3 / 2 incorporate the couplings 6P 3 / 2 <-» (6-40)Si/ 2 and 
6P3/2 (5-42)1)3/2,5/2 ■ The energies of the levels with 
the principal quantum number n < 25 are taken from 
[52I ]. The energies of the levels with the principal quan- 
tum number n > 26 are provided by Arora and Sahoo 
[53I I . The reduced matrix elements for the transitions 
6S1/2 *H> (6-15)P 1 / 2 3 /2 are taken from (54|. The reduced 
matrix elements for the transitions 6P 3 / 2 (6— 10)6'i/ 2 
and 6P 3 / 2 O (5-8)D 3 / 2 5 / 2 are taken from [21[. The re- 
duced matrix elements for transitions to highly excited 
states are provided by Arora and Sahoo [53| . These data 
were calculated by using the relativistic all-order method, 
which includes single and double excitations [21, 54]. The 



calculations for cesium were done in the same way as for 
rubidium [3l|. The full set of parameters we used in our 
numerical calculations is given in Appcndix[DJ The states 
whose energy differences from the ground state are larger 
than the cesium ionization energy of 31406 cm -1 provide 
a discrete representation of the continuum, similar to the 
calculations of Ref. 55] for lithium. We add the contri- 
bution of the core, equal to 15.8 a.u., to the results for the 
scalar polarizabilities [2l|. The polarizabilities are given 
in the atomic unit (a.u.) e^a^/Eh, where ao is the Bohr 
radius and Eh = TO e e 4 /(47reo?i) 2 is the Hartree energy. 

We plot in Fig. [2] the scalar and vector polarizabilities 
a s nJ and ti" njl respectively, of the ground state 65i/ 2 . As 
can be seen, in the region of wavelengths from 400 nm to 
1600 nm, the profiles of both a s nJ and a v nJ have two pairs 
of closely positioned resonances. One pair corresponds to 
the transitions between the ground state 651/2 an d the 
excited state 6Pw 2 [D\ line, wavelength 894 nm) and 
the excited state 6P 3 / 2 (D 2 line, wavelength 852 nm). 
The other pair corresponds to the transitions between 
the ground state 6«S , 1 / 2 and the excited state 7Pi/ 2 (wave- 
length 459 nm) and the excited state 7P 3 / 2 (wavelength 
455 nm) . The effects of the other transitions are not sub- 
stantial in this wavelength region. We note that our nu- 
merical calculations give the values a^ J (6S , 1 / 2 ) — 398.9 
a.u. and a 1 nJ (6S 1 f 2 ) = for the scalar and vector po- 
larizabilities, respectively, of the ground state 615*1/2 of 
atomic cesium in the static limit (u> = 0). The static 
value a s nJ {a>Si/2) — 398.9 a.u. is in agreement with the 
high-precision ab initio theoretical values of 399.8 a.u. 
14 1 and 398.2 a.u. [15| and the experimental value of 
401 a.u. (Hi. 
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FIG. 2. Scalar (a) and vector (b) polarizabilities a s nJ and 
a° nl , respectively, of the ground state 6S 1 j 2 of atomic cesium 
as functions of the light wavelen gth A. The data of this figure 
is provided as electronic files in [33]. 

We plot in Fig. [3] the scalar, vector, and tensor polar- 



izabilities a 



, and a„j, respectively, of the excited 



state 6P 3 / 2 . The figure shows that all the three com- 
ponents have multiple resonances. The most dominant 
resonances are due to the transitions from 6P 3 /2 to (6- 
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8)51/2 and (5-8)133/2, 5/2- We note that our numerical 
calculations give the values a s nJ {QPj,/ 2 ) — 1639.6 a.u., 
<j(6^3/ 2 ) = 0, and alj(QP 3/2 ) -260.4 a.u. for the 
scalar, vector, and tensor polarizabilities, respectively, 
of the excited state 6P3/2 of atomic cesium in the static 
limit (uj = 0). The static values a 8 n j(6P 3 / 2 ) — 1639.6 a.u. 
and a^j(6P3/2) — —260.4 a.u. are in agreement with the 
high-precision ab initio theoretical values of 1650 a.u. 
and —261 a.u. [2l|, respectively, and with the experi- 
mental values of 1641 a.u. and —262 a.u., respectively 



3 

CO 

n, 
o 

« 



N 



40000 



o 




-40000 



400 600 800 1000 1200 1400 1600 
Wavelength of light (nm) 



FIG. 3. Scalar (a), vector (b), and tensor (c) polarizabilities 
a s nJ , oi^j, and a^j, respectively, of the excited state 6P3/2 
of atomic cesium as functions of the light wavelength A. The 
data of this figure is provided as electronic files in [32I ]. 
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FIG. 4. Polarizabilities of the ground state 6S1/2 (red color) 
and the excited state 6P3/2 (blue color) of atomic cesium in 
the region of blue-detuned wavelengths from 680 nm to 690 
nm. The scalar, vector, and tensor components oi a n j, ctnj, 
and qTj are shown by the solid, dashed, and dotted curves, 
respectively. 
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In order to display certain details, we plot in Figs. 0] 
and [5] the polarizabilities a s nJ (solid lines), a v n j (dashed 
lines), and <x[j (dotted lines) of the ground state 6<Si/ 2 
(red color) and the excited state 6P 3 / 2 (blue color) in 
two specific regions of wavelengths. Figures [4] and [5] show 
that the crossings of the scalar polarizabilities a£j (6Si/ 2 ) 
and a^j(6P 3 / 2 ) of the ground and excited states, re- 
spectively, occur at the blue-detuned magic wavelength 
Xb — 686.3 nm [2(| and the red-detuned magic wave- 
length \r ~ 935.2 nm [l9[. Here, red and blue refer to 
the detunings with respect to the P-line transitions. We 
observe from Figs. 0] and [5] that the magnitude of the 
vector polarizability a v nJ is, in general, substantial com- 
pared to that of the scalar polarizability a & nJ - Due to this 
fact, the vector polarizability can contribute significantly 
to the Stark shift when the polarization of the field is not 
linear. 

Comparison between Figs. Ufa) and HJb) shows that 
at the wavelength X v ~ 880.2 nm, which lies between the 
D\ and D 2 lines, the scalar polarizability a s n j(6Si/2) of 



FIG. 5. Same as Fig. [4] but in the region of red-detuned 
wavelengths from 930 nm to 940 nm. 



the ground state is vanishing while the vector polariz- 
ability a^ J (6S , 1 / 2 ) of this state is significant (see Fig. O 
[38l448l |. At this specific wavelength, the ac Stark shifts 
of the sublevels of the atomic ground state are just the 
Zccman-like shifts caused by a fictitious magnetic field 
B . In other words, when specified in the quantization 
coordinate frame {x, y, z} with the axis z parallel to the 
direction of the vector product i[£* x 5], the sublevels 
\FM) of the ground state will be shifted by an amount 
proportional to (— \) F Mi[£* x £]. 

We note that the detunings A2 = Ui v — lod 2 and 
Ai = uj v — ud 1 of the pure-vector-shift (scalar-shift- 
cancellation) frequency u v = 2irc/\ v from the D 2 - and 
Pi-line transition frequencies lud 2 and wd i; respectively, 
are such that A2/A1 = —2.03 ~ —2, in agreement with 
the results of Refs. [1^, |4l|, H2|. In order to understand 
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spectively, of the ground state 6S1/2 of atomic cesium for light 
wavelengths in the region from 870 nm to 890 nm. 



this feature, we make a few additional approximations for 
the scalar polarizability a s nJ in the case where the level 
\nJ) is the ground state nS\/2- We keep only the excited 
levels \n'J') = nP 3 / 2 and \n J') = nP X j 2 in the sum over 
n'J' in Eq. (jTTT) . In the framework of the rotating- wave 
approximation, we neglect the counter-rotating term con- 
taining oJn'j'nj + u in Eq. ([TT]) . We also neglect -fn'j'nj 
in the denominator of the co-rotating term containing 
Wn'J'nJ ~ w - When we insert the result into the first 
expression in Eqs. (|16[) . we obtain the following approx- 
imate expression for the scalar shift of the ground state: 



a n.j = - 



|(nP 3/2 ||d||nS 1/2 )| 2 |(nP 1/2 ||d||nS 1/2 )| 5 



6hA 2 



6?iAi 



It is clear that a s nJ = when 



A, |(nP 3/2 ||d||7iff 1/2 )| 2 
Ar |(nP 1/2 ||d||nS 1/2 )| 2 ' 

With the help of the formula [H] 



(30) 
(31) 



(nJ\\d\\n'J'\ 



\L+S+J'+1 



V(2J+1)(2J' + 1) 



*\ L> S L \(nL\\d\\n'L'), (32) 



we find 



|(nP 3/2 ||d||nff 1/2 ) 
\{nP 1/2 \\d\\nS 1/2 ) 



= 2. 



(33) 



This explains why the relation A 2 ~ — 2Ai is observed 
for the position of A„ in the case of Fig. HI The deviation 
of the ratio A 2 /Ai from the value of —2 is due to several 
reasons. The first reason is that a large number of transi- 
tions are included in our numerical calculations. The sec- 
ond reason is that the counter-rotating terms are taken 
into account in our calculations. The third reason is that 



we used the experimental values |(6P3/ 2 ||d||6S'i/ 2 )| = 
6.324 a.u. and |(6P 1/2 ||d||65 1/2 )| = 4.489 a.u., with 
the ratio |(6P 3/2 ||d||65 1/2 )| 2 /|(6P 1/2 ||d||65 1/2 )| 2 = 1.98. 
The deviation of this ratio from the value of 2 is due to 
relativistic effects 154. 



III. SUMMARY 

We provided a concise, yet comprehensive compilation 
of the general theoretical framework required for calcu- 
lating the polarizability of the states of multilevel atoms 
in light fields with arbitrary polarization. Special em- 
phasis is placed on the interpretation of the vector light 
shift as the result of the action of a fictitious magnetic 
field. We exemplarily applied the presented formalism 
to atomic cesium and calculated the scalar, vector, and 
tensor polarizabilities of the states associated with the 
D2-line transition. Using these results, we highlighted 
points of experimental interest such as the red- and blue- 
detuned magic wavelengths and a wavelength at which 
the scalar light shift of the ground state vanishes while 
the vector light shift is substantial. The underlying set 
of atomic data for the calculations of the polarizability of 
cesium is explicitly given in tabular as well as electronic 
forms. By providing all general tools and definitions in a 
single source and by discussing their respective range of 
validity, our work should facilitate the theoretical model- 
ing of the light-induced potentials experienced by atoms 
in complex far-off-resonance optical fields, encountered, 
e.g., in nonparaxial or near- field optical dipole traps. 
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Appendix A: ac Stark shift of a two-level atom 
interacting with a far-off-resonance light field 

We consider a two-level atom interacting with a far-off- 
resonance light field. Let \a) and \b) be the bare eigen- 
states of the atom, with unperturbed energies E a = hu a 
and Eb = hujf,, respectively, and let w, £, and u be the 
frequency, the complex amplitude, and the complex po- 
larization vector, respectively, of the light field. The elec- 
tric component of the light field is given by Eq. J3|). The 
interaction between the atom and the field is given, in 
the dipole approximation, by Eq. (U). The evolution of 
the off-diagonal density-matrix element pb a of the atom 



9 



is governed by the equation 

Pba = ~i(Ub - Wo - ijba/2)Pba ~ T^T 
+ E*U*e l " t )-d ba (Pbb~ Paa), 



(Sue" 



(Al) 



where dfc a = (6]d|a) is the matrix element of the electric 
dipole operator d = d& |&)(a| + d a b\a)(b\ and ■jba is the 
linewidth of the atomic transition \b) f-> \a). In general, 
we have jba = 7&+7a; where 75 and j a are the decay rates 
of the populations of the levels \b) and \a), respectively. 

We assume that the atom is initially in the level |o), 
which can be, in general, higher or lower than the level 
| b). When the magnitude of the detuning u> — \ui b — ui a \ 
is large compared to the atomic decay rate "fba and to 
the magnitude of the Rabi frequency £1 = d ba £/h, we 
have pbb — and p aa — 1- We use the ansatz pb a — 
Pba e ~ %u}tjr Pba e%U)t an d assume that p\ a and p^ a vary slowly 
in time. Then, we find 



Pba 



Pba 



£u ■ d ba 



1 



2H UJ ba - W - ijba/% 
£*U* ■ d ba 1 



27; 



UJ ba + U - i^ba/2 



(A2) 



The induced dipole is given by p = (d) = db a p a b + 
dabPba- It can be written in the form p = (pe~ luJt + 
p*e lwt )/2, where p = 2(d ba p b ~* + d ab p£ a ) is the envelope 
of the positive frequency component. We find 



£u ■ d ba 



1 



h UJba - OJ - ijba/2 

£u ■ d ab 1 



UJba+UJ + i^ba/2 



(A3) 



The ac Stark shift 5E a of the energy level \a) is the 
time-averaged potential of the induced dipole moment p 
interacting with the driving electric field E and is given 
by 



1- 



1, 



SE a = --p(i) • E(t) = --Re[p ■ £*u*}. (A4) 

Here, the factor of 1/2 accounts for the fact that the 
dipole moment is induced. Inserting Eq. (IA3|) into Eq. 
yields 



SE « = - l -lT Re (— 

4h \uj ba 



|u • d ha | 



|u-d a6 | 2 



■-ijba/2 UJ ba +UJ + ij ba /2 

(A5) 

We emphasize that Eq. (|A5|) is valid for an arbitrary po- 
larization of the light field. When we generalize Eq. (|A5j) 
to the case of a multilevel atom, we obtain Eq. ([5]) . 



Appendix B: ac Stark interaction operator and 
components of the dynamical polarizability 

In this Appendix, we present the details of the deriva- 
tion of the expressions for the ac Stark interaction oper- 
ator V EE and the dynamical scalar, vector, and tensor 



polarizabilities [see Eqs. (|9|)- ([TTj) ]. For this purpose, we 
follow closely Ref. Hf. 

We use the Cartesian coordinate frame {x,y,z}. We 
introduce the notations 



A-i = (A x - iA y )/V2, 
A Q =A Z , 

Ai = -(A x + iAy)/V2 



(Bl) 



for the spherical tensor components of an arbitrary com- 
plex vector A = {A Xl A y , A z }. In terms of the tensor 
components A q = A\ q , with q = —1,0,1, the vector 
A = Ai is an irreducible tensor of rank 1. We intro- 
duce the notation {A ® H}k for the irreducible tensor 
products of rank K = 0, 1, 2 of two arbitrary vectors A 
and B. The q component of the tensor product {A®H}k 
is defined as 

{A ® B}^ = C? q q llq2 A qi B q2 , (B2) 



where 



jimij 2 m 2 



f-i)*-' a+m v/2T+T( 31 32 3 

v 1 mi m 2 — m 



(B3) 

is the notation for the Clebsch-Gordan coefficients. More 
general, an irreducible tensor product of two irreducible 
tensors Uk x and Vk 2 is defined as the irreducible ten- 
sor {Uid ® Vk 2 }k of rank K whose components can be 
expressed in terms of Ux iqi and Vk 2<12 according to 

{U Kl ® Vk 3 }k 9 - £ C^ qiK2q U Kiqi V K2q2 , (B4) 



with K = \K t - K 2 \, \Ki - K 2 \ + 1, K x +K 2 - 1, K 1 + 
K 2 and q = -K, -K + 1, . . . , K - 1, K. Meanwhile, the 
scalar product of two irreducible tensors Uk and Vk is 
defined as 



{u K -v K ) = J2(- 1 ) 9u K, q v K , 



(B5) 



When we use the formula [58j 
(A-B)(A'-B')= (-1) A '{A®AV-{B®B'} K , 



if =0,1, 2 



(B6) 

which is valid for commuting vectors, we can change the 
order of coupling of the operators in Eq. (jSJ) to obtain 



l^l 2 



r"--^ ]T (-i) K {u*®u} K -[{d®n + d} K 

if =0,1, 2 



{-l) K {d®K-d} K }. 



(B7) 



In deriving the above equation we have employed {u (g> 
u*} K — (— l) A '{u* ® u} K - When we use the definition 
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(|B5p for the scalar product of tensors, we can rewrite Eq. 
(lETT) as 



A=0,l,2 



x [{d®72 + d}jr,-« + (-l) A {d®^-d}x,- 9 ]. (B8) 

The explicit expressions for the compound tensor compo- 
nents {u* ® u}^ g , which appear in Eqs. (IB7|) and (|B8I) . 
are 



{u* ® u} ,o 

{u* ® u}i,o = 
{u*(»u}i,i = - 
{u*®u}i,_i = 



V3' 



(B9) 



M 2 -k-il 2 

UqU* + UqU-i 

V2 : 



(BIO) 



and 



{u* <g) u} 2 , = 
{u* (g) u} 2 ,i = - 



{u* ® u} 2 



3|u | 2 -l 



V2 ' 

{U* <g> U} 2 , 2 = -U\U*_ x , 

{u* (gu} 2 ,- 2 = -u-ml. (Bll) 

The operators 7?.+ and 1Z- in Eqs. (|B7|) and (|B8|) are 
given by Eqs. ©. In our treatment given below, the ba- 
sis states | a) and \b) in Eqs. ([7J are taken from the F ba- 
sis states \nJFM), with unperturbed energies U! n JFM = 
uj n j and spontaneous decay rates j n JFM = 7nj- 

Let V^ F , M , ee ((nJ)FM|V r££; |(nJ)F'M') be the 
matrix elements of the Stark interaction operator V EE 
in the atomic hfs basis {\(nJ)FM)} for a fixed set of 
quantum numbers nj. From Eq. (IB8|) . we find 



9 

(B12) 



A=0,l,2 



A ^FMF'M'i 



where 



-9 



FMF'M 1 ~ / j °lai 1 
9192 

x E (nJFM|d gi |n"J"F"M") 



n" J" F" M" 



x (n W'M"K 2 \nJF'M'}n^ JltnJ , (B13) 



with 



/'' <C - ^„., - -c- - /-„»./»„.//2 



(-1) 



A' 



~ w nJ + w + «7„» j»„j/2 



(B14) 



According to the Wigner-Eckart theorem [53, the de- 
pendence of the matrix elements (nJFM\TK q \n' J' F 1 M') 
of tensor component operators T Xq on the quantum num- 
bers M, M', and q is entirely included in the Wigner 3-j 
symbol, namely, 



{nJFM\T Kq \n'J'F'M' x 



(-1) 



F — M 



F K F' 



M a M> )(nJF\\T x \\n'J , F'). (B15) 



Here, the invariant factor 



, F -ul F K F< 
-M q M' 



(nJF\\T K \\n'J>F>) = £ (-1 

MM'q 

x (nJFM\T Kq \n'J'F'M'} (B16) 



is the reduced matrix element for the set of tensor com- 
ponent operators Ta<j, with the normalization convention 

\(nJF\\T K \\n'J'F')\ 2 = ^ \(nJFM\T Kq \n'J'F'M')\ 2 

MM'q 

(B17) 

and the complex conjugate relation 

(nJF\\T K \\n'J'F')* = (-1)^' (n'J'F'\\T K \\nJF). 

(B18) 

Since the electric dipole d is a tensor of rank 1, the ap- 
plication of the Wigner-Eckart theorem to the matrix 
elements (nJFM\d q \n' J'F'M') of the spherical-tensor- 
component operators d q of the electric dipole gives 



{nJFM\d q \n' J'F'M') 



= (-1) 



1 

F-M 



F 1 F' 

-M q M' 



{nJF\\d\\n'j'F'). (B19) 



The invariant factor (nJ F\\d\\n' J' F') is the reduced ma- 
trix clement for the electric dipole operator d. With the 
help of Eq. (|B19|> . we can rewrite Eq. (|B13|) as 



O 



Kq 

FMF'M' 



E (nJF\\d\\n"J"F") 



n"J"F" 



x (n''J''FH|d||nJFO^„ rlJ ^ F F ; M ,, (B20) 

where 

M^;, M , = V2KTT E hl) F + F»-M-M»- q 
9192M" 

1 1 K \f F IF" 
qi qi q )\ -M qi M" 

F" 1 F' 



-M" q 2 M' 



(B21) 
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When we use the symmetry properties of the 3-j symbol 
and the sum rule [581 



X 



+J5+J6— ITU— ms— me ( 35 Jl 36 
7775 — nil — 777-6 



J6 J2 J4 
777g —7772 — 7774 



J4 J3 J5 
777,4 — "73 — fTT-5 



Jl 


h 


h \ 




n 


33 1 


mi 


77l 2 


m 3 J 


{Ji 


k 


ie J 



we find 

JV FMF'M< 



(-1) K+F+M V2K + T 



1 X 1 

F F" F' 



, (B22) 



F A F' 
-M -q M' 

(B23) 



We now insert Eq. (|B23|) into Eq. (|B20|) and then insert 
the result into Eq. (|B12[) . Then, we obtain [H[ 



EE 

FMF'M- 



Kg 



if =0,1, 2 



F K F' 



where 

^ = (-1) A ™'V2^TT £ { 

n"J"F" 

X (77jf||d||77"j"F")(77"j"^"||d||r7j^)<». ) / »„. / 

(B25) 

are the reduced scalar (A = 0), vector (if = 1), and 
tensor (A = 2) polarizability coefficients for the hfs levels 
within a fine-structure manifold n J. 

For the tensor component operators Tx q that 
do not act on the nuclear spin degrees of free- 
dom, the dependence of the reduced matrix element 
(nJIF\\T K \\ri J'I'F') on F, F', I, and F may be fac- 
tored out as p58j 

(nJIF\\T K \\ri J'I'F') = 5 I r(-l) J+I+F ' +K 

x V(2F+1)(2F' + 1){ j, * F j }<rcJ||TV||rc'j'}. 

(B26) 

Since the electric dipole d of the atom does not couple 
to the nuclear degrees of freedom and is a tensor of rank 
1 , the use of Eq. (|B26|) for the case Tk = d yields 



1 K 1 

F F" F' 



iJF\\d\\n'J'F 



' J' T?'\ — 



iJIF\\d\\n'J'IF') 



= {-l) J+I+F ' +1 ^(2F + 1)(2F' + 1) 

F IF' 
J' I J 



iJ\\d\\n'J r 



(B27) 



Substituting Eq. (|B27)) into Eq. (|B25|) yields 

S K ) -I y\K+J+2I+F+2F' 



<7ff' = (-l) K+J+ ' zl+m W(2F+l)(2F> + l) 

X (-l) J "(«J||d||77"j")(77"j"||d||77j)<^, nJ 



n" J 

X 



V2A + 1^(-1) F "(2A" + E 



F 1 K 





1 


F" \ 




i 


A" \ 




I 


J J 




i 


J J 



(B28) 
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J5 J 


1 ie 




J8 J 



The summation over A" in Eq. (IB28I) can be performed 
using the formula [HI 



]>>im+i) 



J7 

ia ji 39 



The result is 

S K ) - 



(—\\ ~h ~i 2 _ & _ -? 7 _ -? 9 

(B29) 



1* 


is 


J9 1 


f J5 


is 


i'g 1 


i 36 




h J 


\ J> 


J2 


i3 J 



Q 



n.JFF' 
X 



= (-1) /+F '-V(2F+1)(2F' + 1) 
53<nJ||d||n"J")<n"J"||d||nJ)72W„ Bj 



n"J" 

x V2A + 1 





A' 






A" 







./" 






I 





(B30) 



When we insert the explicit expression (|B14[) into the 
above equation, we get [24| 



x nJFF 



, =(-l) J+I+F ' +K v / {2F + 1)(2A' + 1) 



F K F< i (K) 
J I J ' " 



(B31) 



where 



(if) 



= (-1) 2J+k+1 n /2ATiW 



1 A 1 

J J 1 J 



X (77j||d||77'j')(77'j'||d||nJ) 

1 / 1 

x -Re — 

h \LOn'J' - UnJ -(*>- Hn'J'nJ/2 



(B32) 



(-1) 



K 



OJn'J' - VnJ + ^ + 77„',/'„j/2 



are the reduced scalar (A = 0), vector (A = 1), and 
tensor (A = 2) polarizabilities for the Stark shift of 
the fine-structure level \nJ). When we substitute Eq. 
(|B3T|) into Eq. ([B24| . we obtain Eq. (fTU]). Since 
(nJ\\d\\n'J') = (-l) J - J '(77'J'||d||77j), Eq. (jB32l) can be 
rewritten as Eq. (jll[) . The explicit expressions for the 
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reduced polarizabilities o&fp are 
2 



x nJ 



==£Kn'J'||d||nJ)| s 

> n'J' 



a 



(i) 

nj 



v jWjjw - ^ 2 + 7n' J'nj/ 4 ) 

2X "^' ( j ), ))|(n'J'||d||nJ)| 2 



; E( _ 1) . + .r| 

n'J' ^ 



^I'J'nJ - - ll'J'nj/^) 



(2) 
nJ 



("I'J'nJ -^ 2 + llmjl^) 2 + ll'J'n.^ 2 ' 

' -'- A ' 2 U|(n'J'||d||nJ)| 2 



^ 2 +l 2 n'J'nj/^ 2 +^'J'nJ^ 



n' J'nJ 



(B33) 



When we neglect the linewidths jn'j'nj, Eqs. (|B32[) and 
()B33j) come to full agreement with the results of Ref . [24| . 

We emphasize that, in the above calculations for the 
polarizabilities, we used the approximation io n jFM — 
uj n j, that is, we neglected the effect of the hfs splitting 
on the polarizabilities. This approximation is consistent 
with the perturbation theory scheme used in our case 
where the hfs splitting and the ac Stark shift are consid- 
ered to be small perturbations of the same order. If the 
hfs splitting is much larger than the ac Stark shift, we 
can consider only the ac Stark shift as a perturbation. In 
this case, the operator for the ac Stark shifts of sublevels 
of a hfs level \nJF) is given by expression 24 [ 



V EE = -~\E\ 



1 



+ a 



, [u*xu]-F 
4'-' V F 2F 
T 3[(u* • F)(u • F) + (u • F)(u* • F) 



2F 2 



2F(2F - 1) 



(B34) 



with the matrix elements 



V, 



EE 
MM' 



\£\ 2 



Kq 



K=0,l,2 

X(-1) F " M (_M K q M>H K) > ( B35 ) 



F K F 



where 



„(°) 



V3(2F+1) 



2F 



(F+1)(2F + 1) 



(1) 
F ' 



2F(2F- 1) 



a ' F = -^3(F + l)(2F+l)(2F + 3) a F> (B36) 



and 



n' J' 



af ] = {-1) K+F+1 (2F+1)V2KTTJ2 \(^'J'M\nJ)\' 

xE(-l) F ' (2^ + 1)1 
pi v 





A' 






l 






F' 











1 / 1 

x -Re ; j- 

n \U} n 'J'F' — U n JF — OJ — l-fn'J'F'n.JF/2 



(-1) 



A" 



^n'J'F' — U n jp + LO + ijn'J'F'n.JF/2 



(B37) 



We note that, in the framework of the validity of Eqs. 
(jB34|) and (|B37)) . the different hfs levels F = I + 1/2 
and F — I — 1/2 of the ground state have different scalar 
polarizabilities. This means that, when the hfs splitting 
is taken into account in the expression for the atomic 
transition frequency u) n iji pi—uj n jp, a nonzero differential 
scalar Stark shift between the hfs levels of the ground 
state may occur. 



Appendix C: Additional magnetic field 

We consider the presence of a weak external real mag- 
netic field B. The Hamiltonian for the interaction be- 
tween the magnetic field and the atom is [l| 



V B = fi B g n j{3-B). 



(CI) 



It can be shown that the matrix elements of the operator 
V B in the basis {\FM)} are given by the expression 



VP 



FMF'M' 



VBgnj(-l) J+I - M VJ(J+ 1)(2J+ 1) 

x v/(2F+l)(2F' + l)/ 1 ' 1 



^ ("Wl M I -M> ) lC - >) 

9=0, ±1 



J I J 

F' 
q —M' 



Here, 5_ x = (B x - iB y )/y/2, B = B z , and B 1 = 
— (B x + iB y )/^/2 are the spherical tensor components of 
the magnetic induction vector B = {B x , B y , B z }. We 
note that Eq. (|C2[) is valid for an arbitrary quantization 
axis z. 

When F is a good quantum number, the interaction 
operator (|Clj) can be replaced by the operator 



V 1 



Mb5«jf(F • B) 



(C3) 



In the absence of the light field, the energies of the 
Zeeman sublevels are huj n jpM = Tiu n jp + [iBQnjF 

BM. 

Here, huj n jp is the energy of the hfs level \nJF) in the 
absence of the magnetic field and M — — F, . . . , F is the 
magnetic quantum number. This integer number is an 
eigenvalue corresponding to the eigenstate \FM)b of the 
projection F ZB of F onto the zb axis. In general, the 
quantization axis z may be different from the magnetic 
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field axis zb and, consequently, \FM) may be different 
from \FM)b- In order to find the level energy shifts, 
we must add the magnetic interaction operator V B to 
the combined hfs-plus-Stark interaction operator (J5| and 
then diagonalize the resulting operator. 

Appendix D: Atomic level energies and reduced 
matrix elements 

The full set of parameters used in our numerical calcu- 
lations is shown in Tables HHVl The data is also provided 



as electronic files in (32|. The reduced matrix elements 
are given in the atomic unit eao of the electric dipole. The 
energies of the levels with the principal quantum number 
n < 25 are taken from (52[. The energies of the levels 
with the principal quantum number n > 26 are provided 
by Arora and Sahoo [53| . The reduced matrix elements 
for the transitions 6S1/2 O (6-15)Pi/2.3/2 are taken from 
[54j ]. The reduced matrix elements for the transitions 
6P3/2 (6-10)5i/ 2 and6P 3 / 2 <h- (5-8)£> 3 / 2 ,5/2 are taken 
from [2lj . The reduced matrix elements for transitions 
to highly excited states are provided by Arora and Sahoo 
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TABLE I. Energies of levels nPi/2 and reduced matrix ele- 
ments for transitions between levels nP 1 / 2 and 65V/2 m atomic 
cesium. The data of this table is provided as an electronic file 
in [12. 



nPi/2 level 


Energy 


\(nP 1/2 \\d\\6S 1/2 )\ 




(cm" 1 ) 


(a.u.) 


6P1/2 


11178.27 


4.489 


7Pl/2 


21765.35 


0.276 


8P1/2 


25708.84 


0.081 


9Pl/2 


27637.00 


0.043 


10P 1/2 


28726.81 


0.047 


HPl/2 


29403.42 


0.034 


12P1/2 


29852.43 


0.026 


13Pl/2 


30165.67 


0.021 


WP1/2 


30392.87 


0.017 


15Pl/2 


30562.91 


0.015 


16P1/2 


30693.47 


0.022 


17Pl/2 


30795.91 


0.023 


18P1/2 


30877.75 


0.019 


19P 1/2 


30944.17 


0.010 


20P 1/2 


30998.79 


0.035 


2IP1/2 


31044.31 


0.002 


22P 1/2 


31082.60 


0.027 


23P 1/2 


31115.12 


0.002 


24P 1/2 


31142.97 


0.000 


25P 1/2 


31167.02 


0.031 


26P 1/2 


79752.85 


0.045 


27P 1/2 


123985.9 


0.044 


28P 1/2 


201030.4 


0.040 


29P 1/2 


334252.0 


0.035 


30P 1/2 


563660.5 


0.029 


31Pi/2 


958103.9 


0.022 


32P 1/2 


1636312 


0.015 


33P 1/2 


2804845 


0.010 


34P 1/2 


4827435 


0.007 


35P 1/2 


8334762 


0.004 


36P 1/2 


14384903 


0.002 


37P 1/2 


19069034 


0.000 


38P 1/2 


24735902 


0.001 


39P 1/2 


42263044 


0.001 


40P 1/2 


71483200 


0.000 
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TABLE II. Energies of levels nP 3 / 2 and reduced matrix ele- 
ments for transitions between levels nP 3 / 2 and 65 , 1 / 2 m atomic 
cesium. The data of this table is provided as an electronic file 
in [12. 



nP 3 / 2 level 


Energy 


KnP 3/2 ||d|l6S 1/2 }| 




(cm- 1 ) 


(a.u.) 


6P3/2 


11732.31 


6.324 


7P3/2 


21946.40 


0.586 


8P3/2 


25791.51 


0.218 


9P 3 /2 


27681.68 


0.127 


IOP3/2 


28753.68 


0.114 


HP3/2 


29420.82 


0.085 


I2P3/2 


29864.34 


0.067 


13P 3 /2 


30174.18 


0.055 


WP3/2 


30399.16 


0.046 


15P 3 /2 


30567.69 


0.039 


I6P3/2 


30697.19 


0.062 


17P 3 /2 


30798.85 


0.065 


I8P3/2 


30880.12 


0.062 


19P 3 /2 


30946.11 


0.039 


20P 3/2 


31000.40 


0.112 


2IP3/2 


31045.66 


0.014 


22P 3/2 


31083.77 


0.119 


23P 3/2 


31116.09 


0.081 


24P 3/2 


31143.84 


0.004 


25P 3/2 


31167.74 


0.032 


26P 3/2 


82736.05 


0.001 


27P 3/2 


129247.4 


0.017 


28P 3/2 


210254.6 


0.025 


29P 3/2 


350320.8 


0.029 


30P 3/2 


591569.5 


0.029 


3IP3/2 


1006618 


0.023 


32P 3/2 


1720992 


0.017 


33P 3/2 


2953614 


0.012 


34P 3/2 


5089655 


0.008 


35P 3/2 


8794667 


0.005 


36P 3/2 


15185971 


0.003 


37P 3/2 


19069122 


0.000 


38P 3/2 


26120330 


0.002 


39P 3/2 


44625560 


0.001 


40P 3/2 


75454848 


0.001 
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TABLE III. Energies of levels nSi/ 2 and reduced matrix ele- 
ments for transitions between levels nSi/ 2 and 6P3/2 m atomic 
cesium. The data of this table is provided as an electronic file 
in [13. 



nSi/2 level 


Energy 


|<nS 1/2 ||d|l6P 3 /2>| 




(cm- 1 ) 


(a.u.) 


6S1/2 





6.324 


7<S'i/2 


18535.53 


6.470 


8<5'i/2 


24317.15 


1.461 


9S 1/2 


26910.66 


0.770 


IOS1/2 


28300.23 


0.509 


HSi/a 


29131.73 


0.381 


I2S1/2 


29668.80 


0.297 


13S'i/2 


30035.79 


0.241 


us 1/2 


30297.64 


0.219 


1551/2 


30491.02 


0.234 


I6S1/2 


30637.88 


0.251 


17Si /2 


30752.03 


0.259 


I8S1/2 


30842.52 


0.239 


19Si/ 2 


30915.45 


0.376 


2OS1/2 


30975.10 


0.213 


2IS1/2 


31024.50 


0.349 


22Si/2 


31065.88 


0.472 


23Si/2 


31100.88 


0.578 


24Si/ 2 


31130.75 


0.026 


25Si/ 2 


31156.44 


0.464 


26S1/2 


74266.21 


0.325 


27Si/ 2 


114033.8 


0.211 


28S1/2 


184506.6 


0.132 


29Si/2 


308495.8 


0.081 


30Si/2 


525668.9 


0.048 


3IS1/2 


905144.4 


0.027 


32Si/2 


1567994 


0.015 


33Si/2 


2728852 


0.008 


34Si/ 2 


4773681 


0.004 


35Si/2 


8387611 


0.002 


36S'i/2 


14752595 


0.001 


3751/2 


19071256 


0.000 


3851/2 


25891848 


0.001 


39Si/ 2 


45209868 


0.000 


40Si/2 


78196928 


0.000 
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TABLE IV. Energies of levels nD 3 / 2 and reduced matrix 
elements for transitions between levels nD 3 / 2 an d 6P 3 / 2 m 
atomic cesium. The data of this table is provided as an elec- 
tronic file in [32l ]. 



nD 3 / 2 level 


Energy 


|(nL»3 /2 ||d||6P 3 /2>| 




(cm- 1 ) 


(a.u.) 


5£>3/2 


14499.26 


3.166 


6-D3/2 


22588.82 


2.100 


7D 3 /2 


26047.83 


0.976 


8D3/2 


27811.24 


0.607 


9D 3 /2 


28828.68 


0.391 


IOD3/2 


29468.29 


0.304 


H-D3/2 


29896.34 


0.246 


I2D3/2 


30196.80 


0.211 


13D 3 /2 


30415.75 


0.215 


14L» 3/2 


30580.23 


0.234 


15L» 3 /2 


30706.90 


0.248 


I6D3/2 


30806.53 


0.256 


17i?3/2 


30886.30 


0.269 


18D3/2 


30951.15 


0.204 


19D3/2 


31004.59 


0.397 


20D3/2 


31049.14 


0.012 


21D3/2 


31086.68 


0.482 


22L» 3 /2 


31118.60 


0.021 


23D3/2 


31145.97 


0.491 


24L» 3 /2 


31169.61 


0.438 


25L» 3 /2 


31190.18 


0.344 


26D 3 /2 


53740.33 


0.028 


27L» 3 /2 


67538.40 


0.244 


28L» 3 /2 


89286.16 


0.158 


29L» 3 /2 


122155.2 


0.095 


3OD3/2 


171666.2 


0.053 


3ID3/2 


246044.2 


0.027 


32L» 3 /2 


357540.4 


0.011 


33D 3 /2 


524435.1 


0.002 


34L» 3 /2 


774031.9 


0.002 


35L» 3 /2 


1147062 


0.004 


36D 3/2 


1704131 


0.004 


37L» 3 /2 


2535326 


0.004 


38L> 3 /2 


3774912 


0.003 


39L» 3 /2 


5623814 


0.002 


4OD3/2 


8382140 


0.002 


4ID3/2 


12493172 


0.001 


42D 3 /2 


18603884 


0.001 
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TABLE V. Energies of levels nD 5 / 2 and reduced matrix el- 
ements for transitions between levels nD 5 / 2 and 6P 3/ / 2 m 
atomic cesium. The data of this table is provided as an elec- 
tronic file in [32l ]. 



nD 5 /2 level 


Energy 


|<nD B /a||d||6Pa/a)| 




(cm- 1 ) 


(a.u.) 


5D 5 /2 


14596.84 


9.590 


6D5/2 


22631.69 


6.150 


7£>5/2 


26068.77 


2.890 


8D5/2 


27822.88 


1.810 


9D 5 /2 


28835.79 


1.169 


IOD5/2 


29472.94 


0.909 


llAs/a 


29899.55 


0.735 


I2D5/2 


30199.10 


0.630 


13D 5 /2 


30417.46 


0.642 


14D 5 /2 


30581.53 


0.699 


15D 5 /2 


30707.91 


0.741 


I6D5/2 


30807.33 


0.766 


17D 5 /2 


30886.94 


0.798 


I8D5/2 


30951.68 


0.745 


19L»5/2 


31005.03 


0.903 


20D 5/2 


31049.52 


0.840 


2ID5/2 


31087.00 


1.438 


22L» 5 /2 


31118.87 


0.130 


23D 5/2 


31146.20 


1.456 


24D 5/2 


31169.81 


1.288 


25L» 5 /2 


31190.35 


0.149 


26D 5/2 


54508.05 


0.998 


27D 5/2 


69450.32 


0.713 


28L» 5 /2 


92156.01 


0.462 


29D 5/2 


126462.2 


0.280 


30D 5/2 


178123.1 


0.158 


3ID5/2 


255714.2 


0.081 


32L> 5/2 


372021.2 


0.035 


33L» 5/2 


546142.6 


0.008 


34L» 5 /2 


806627.1 


0.005 


35L» 5/2 


1196092 


0.010 


36L» 5/2 


1778024 


0.011 


37L» 5 /2 


2646960 


0.010 


38L» 5/2 


3944018 


0.008 


39L» 5/2 


5880432 


0.006 


40D 5/2 


8771401 


0.005 


4ID5/2 


13082459 


0.003 


42D 5 /2 


19494500 


0.002 



